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ABSTRACT 

Many planets are observed in stellar binary systems, and their frequency may be comparable to 
that of planetary systems around single stars. Binary stellar evolution in such systems influences the 
dynamical evolution of the resident planets. Here we study the evolution of a single planet orbiting one 
star in an evolving binary system. We find that stellar evolution can trigger dynamical instabilities 
that drive planets into chaotic orbits. This instability leads to planet-star collisions, exchange of 
the planet between the binary stars ( "star- hoppers" ) , and ejection of the planet from the system. 
The means by which planets can be recaptured is similar to the pull-down capture mechanism for 
irregular solar system satellites. Because planets often suffer close encounters with the primary on the 
asymptotic giant branch, captures during a collision with the stellar envelope are also possible. Such 
capture could populate the habitable zone around white dwarfs. 



1. INTRODUCTION 

Most of the (> 2000) planets and planet ca ndidates 
discovered to date are hosted by single stars (Howard 
et al. 2011). This bias towards single stars may be pri 



manly due to the avoidance of multiple systems in target 
selection: the limited number of planet discoveries in bi 
nary systems suggests that the planet frequency in mul 
tiples may be comparable to that of single stars |E ggen 



berger et al.||2004t|Bonavita fc D csidcra 2007; Mugrauer 

Neuh auscr 2 009)). ' 
terns (Raghavan et al.||20I0p, a large fraction of all exo 



& Neuhauscr 2009) . Given the frequency of binary sys 



planetary systems may reside in binaries. Here we focus 
on the dynamics of planets in evolving binary systems 
near the peak of the binary separation distribution: from 
- 50 - 100 AU. 
Though various studies have explored the possible for- 



systems (Holman & Wiegert||1999; David et al. 2003 


Mudryk & Wu|200ti |Perets 2010 


Moeckelfe Veras|2012|), 



been thoroughly explored. The influence of stellar evo- 
lut ion on planets in singleton s ystems has been s tudie d 
bylDebes fc Sigurdssonfi J20021 ; jVillaver fc Livio] poM) ) ; 
|Veras et ai] ( |201l| ). |Perets| ( |2010[ ) described qualitatively 
the evolution of planets m evolved binary systems, in- 



cluding both circu mstellar and circumbinary planets. Ve 



ras & Tout ( 2012 ) have recently conducted a quantitative 
study of circumbinary planetary dynamics in evolving 
systems. 

In this paper we focus on the effects of binary stel- 
lar evolution on planetary systems in "satellite" orbits, 
those around one star in the system. We consider the 
influence of slow mass loss via stellar winds, and demon- 
strate that these systems have rich dynamical behavior, 
including planetary ejections, "star hopping" (exchanges 
of planets between the binary stars) , capture of planets 
into stable orbits, collisions with either of the stars, and 
interactions potentially leading to capture into close or- 
bits. A similar mechanism that operates in triple stellar 
systems, dubbed the Triple Evolution Dynami cal Insta- 
bility (TEDI) has been recently explored by Perets & 
Kratterl (120121). 



In the following we explore the outcomes of this insta- 
bility using few-body simulations, and provide a quali- 
tative and semi-quantitative understanding of the prop- 
erties of such systems. We also provide an estimate for 
the fraction of systems that are subject to the planetary 
version of the TEDI, given various assumptions for the 
population of planetary systems in binaries. 

We begin in Section [2] by describing the orbital evolu- 
tion and stability criterion for planets in evolving bina- 
ries. We next explain the capture mechanism for planets 
around the secondary star in the context of the circular 
restricted three body problem (Section [3]). In Section [4] 
we describe a series of n-body experiments that we use 
to explore the outcomes of the chaotic planetary orbits 
induced by stellar mass loss. In Section [5] we present 
probabilities for collision and capture for planets in our 
numerical experiments, and provide a summary of ex- 
pected outcomes. In Section [6] we explore the frequency 
of unstable planetary binary systems. Finally, in Section 
[7] we conclude with several observational predictions for 
next generation surveys. 

2. DYNAMICAL EVOLUTION OF CIRCUMSTELLAR 
PLANETS IN EVOLVING BINARY SYSTEMS 

In order to remain dynamically stable, a triple system 
of gravitating objects requires a hierarchical configura- 
tion with orbits that arc well separated. For satellite type 
planetary-binary systems, there is a critical ratio between 
the semi-major axes of the binary and the planet beyond 
which the system is unstable. The stability boundary 
is quite complex i n detail, set by the overlap of mean 
motion resonances (Wisdom 1980 Mudryk fc Wu|2006| , 
however it can be we ll tit as a function of bin ary mass 
ratio and eccentricity (Holman & Wiegert 1999 ). For cir- 



cular binary orbits, which we consider here, the radius of 
the stable region around the primary is: 



Rn = 0.464-0.38- 



m 2 
mi+m 2 



(1) 



where a* is the stellar semi-major axis, and m\ and mi 
are the masses of the primary and secondary, respec- 
tively. Note that this is identical to Equation (1) of |Hol-| 
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man & Wiegert (1999), but without the quadratic de- 
pendence on eccentricity. These authors find that plan- 
ets with semi-major axis a p > R p are destabilized within 
10 4 binary orb ital periods, (see als o the results of David 
et al.|2 003 and Fat uzzo et al. 12006} who calculate typical 
instability timescaies). The size of the stability region 
around the secondary, R s , is obtained by the substitu- 
tion mi -f-> ni2- Note that the initial eccentricity of the 
planet, which we set t o zero here, is unimp ortant for sta- 



bility considerations (iMudryk fc Wu|2006| ) 

We now consider the impact ot slow (adiabatic) and 
isotropic mass loss from the primary star in the binary, 
where the mass loss timescale is long compared to the 
planetary and stellar orbital periods. Mass loss is the 
byproduct of stellar evolution and occurs through stellar 
winds. Adiabatic mass loss drives the orbits to larger 
semi-major axes: 

a i = W f ai ' (2) 



(Hadjidemctriou 19631, where Mf is the final mass of 
the system alter its evolution, and Mj and a,i are the 
initial mass and initial semi-major axis of the system, 
respectively. This relation derives from a perturbation 
analysis that shows that the the specific angular mo- 
mentum GMa(l — e 2 ) (per reduced unit mass) is con- 
served, and the assumption that mass loss is slow and 
isotropic, which guarantees that eccentricity is also con- 
served. These assumptions are only accurate for planets 
far from their host star (i.e. beyond the tidal effects of 
interaction with the expanding envelope of the evolving 
primary), and where mass transfer is negligible. The 
outcomes of the latter possibility have been discussed 



(e.g. Livio & Sokerl 


1984 


Villaver & Livio 


2007 


Nord- 


haus et al. 2010 


and references therein). 



the planet expand due to mass loss, the change to the 
planet's orbit is much larger than the change in the bi- 
nary's orbit because the mass loss within the planet-star 
system comprises a greater fraction of the enclosed mass. 
Allowing the planet mass to go to zero, the ratio between 
the semi-major axis of the planet to stellar companion 
therefore increases by a factor of: 



a pJ 
a *J 



I 



mi,i 



mij + m 2 



mi. 



m 2 



(3) 



where the subscripts correspond to the system compo- 
nents (1,2, and p for the evolving star, its stellar com- 
panion and its planetary companion, respectively), and i 
and / corresponding to the initial and the final system, 
respectively. Relative orbital expansion due to mass loss 
can shift planets across the stability threshold. 

Combining Equations |T]) and ([3]), we can define the 
approximate stability boundary as a function of the evo- 
lutionary state (i.e. current primary mass) of a given 
binary system. In Figure [lj we show the theoretical evo- 
lution of several planetary orbits as a function of their 
separation ratio for the evolution of a 2M star in a bi- 
nary orbit with a 1M Q companion. In this system, the 
ratio of the planetary to stellar orbit, a p /a*, expands by 
nearly a factor of two as the primary evolves to become a 
white dwarf (hereafter WD). At the same time, the criti- 
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Fig. 1. — Top panel: Evolution of the ratio of planet to stellar 
semi-major axis (thin colored lines) and critical stability bound- 
ary (thick dashed line) for a mass losing 2Mq with solar mass 
companion. Initial semi-major axis ratios which cross the stability 
threshold will enter a chaotic phase of evolution. We show only 
the ~ 2Myr phase where mass loss is most rapid. The stability 
threshold declines more sharply than the mass because the stability 
boundary also shrinks due to the change in fi = m2/(rri2 +mi). 
Bottom panel: the corresponding stellar mass and stellar radius 
evolution for the evolving primary. 

cal boundary for stability around the primary shrinks as 
the primary becomes less massive than its companion. 
We show the evolving mass and radius of the primary as 
well: the stellar model is calculated from the Single Stel- 



lar Evolution (SSE) code of |Hurley etaL ( 2000[ ), which 
evolves the primary from mi = 2.0 —> O.55M . 

Planets with semi-major axis ratios that exceed the 
critical point in Equation ([!]) enter into a phase of chaotic 
evolution. Note that we use chaotic in the technical 
sense, wherein planets with nearly identical initial con- 
diti ons undergo dramatica lly different orbital evolution 
See Mudryk & Wu ( |2006[) for a discussion of measured 
Lyapunov exponents tor the problem without mass loss. 

After becoming unstable, planets can suffer close en- 
counters with either component of the binary. Intrigu- 
ingly, the instability can also lead to the long-term cap- 
ture of planets around the secondary star. We now dis- 
cuss the capture mechanism in more detail. 

3. CAPTURE AND THE JACOBI CONSTANT 

Before describing our numerical results, we begin with 
a brief description of the capture process, in which a 
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planet originally orbiting the primary star is destabilized 
due to stellar mass loss, and then recaptured into a stable 
orbit around the companion. We do so by considering 
the circular restricted three-body problem: coplanar test 
particles moving about two massive bodies on a circular 
orbit. 

Capture occurs through a process reminiscent of the 
suggested pull down m echanism for the capture of 
the Moon (Ruskol 19611 and irregular s atellite capture 
around the gas giants m the Solar Sys tem ( jHeppenheimer 
fc Porco|1977| |Astakhov et al .||2003[). In these cases, the 



body that captures a satellite does so via gaining mass. 
We show here that mass loss by the primary can also 
lead to capture. 

Capture is controlled by the evolution of a particle's 
Jacobi constant, C,j. Recall that the Jacobi constant 
is the only conserved integral of motion in the circular 
restricted three body problem (hereafter CR3BP): 



Cj 



\n r 2 J 



where fix and ^2 ar e the mass ratios of the stars, fj and 
T2 are the separations of the particle from each star, n is 
the mean motion, and the cartesian components are the 
positions and vel ocities of the particle in th e non-rotating 
reference frame (Murray & Dermott 2000). A particle's 
Cj defines the parameter space in which it may orbit the 
binary. The boundaries of the permitted region can be 
visualized through the use of zero- velocity curves. These 
contours are obtained by setting the planet velocities to 
zero and solving Equation Q for the positions in the 
non-rotating reference frame. 

By adding mass loss into the CR3BP, Cj is no 
longer conserved, because /i evolves in time. However, 
Luk'yanov (2009) has shown that this so-called quasi- 
Jacobi integral is mathematically well behaved in the 
case of mass transfer in a binary system, suggesting that 
it still specifies the bounds of the planet's orbit at any 
instant. 

The evolution of Cj due to mass loss dictates whether a 
given particle can pass between the stars, or must remain 
in orbit about one or the other. Figure [2] shows the zero 
velocity curves for a planet at four epochs in the evolution 
of a system with mi = 2.0 — > 0.55M Q , m-i = l.OiW©, 
and initial semi-major axis, a*,i = 90 AU. Figure [3] is 
complementary, and shows the evolution of the planetary 
and stellar separations, as well as the particle's Cj as a 
function of time. 

The fall and rise of Cj shown in both Figures [2] andjs] is 
characteristic of host switching. As shown in Figure [3l we 
can estimate the evolution of Cj analytically for a par- 
ticle's orbit that expands with mass loss as dictated by 
Equation ([5]). At early times, we keep the planet fixed at 
its Keplerian rotation velocity. Once the planet nears the 
instability threshold, the orbit is highly non-Keplerian. 
In this regime we interpolate between Cj calculated for 
two assumptions: a) that the planet remains on a Ke- 
plerian orbit about the primary, and b) that the planet 
switches to a Keplerian orbit at the same separation, 
but about the secondary. The disagreement between the 
analytic estimate and the actual particle value is worst 
where orbits depart most from Keplerian (i.e. in the 
planet hopping region). The relatively better agreement 



at late times suggests that the particle is well described 
by a Keplerian orbit, and that the orbital energy is con- 
sistent with that of the initial orbit. 

In the case shown here (the simulations are described 
in the following section), the planet's Cj at first pro- 
hibits transfer between the two stars, then permits it, 
and then once again forbids it. If a planet is in the sta- 
ble region around the secondary when the curves close 
off, the planet becomes trapped. One can distinguish 
between open and closed zero-velocity curves by com- 
paring the particles' Jacobi constants to the value at the 
first Lagrange point ( £1), or: 



Cj(Lx) 



. 3(4/3)^(2/3) 



1(W3 



(5) 



( |Murray fc Dermott||2000[ ). 

IN ot all systems undergo this complete cycle. For some, 
when mass loss stops, the zero velocity curves remain 
open, and capture is less likely. Note that the stability 
threshold defined in Equation (fil is crossed only after the 
curves have opened: planets can remain on stable orbits 
in the face of open zero-velocity curves. 

We show in Section [5] that the relative value of the the 
average Jacobi_constant of the surviving particles and at 
the Li point, Cj-Cj(Li), correlates with the probability 
of planets being captured. 

4. FEW-BODY SIMULATIONS 

We use few-body simulations to explore the evolution 
of planets in evolving binary systems. We focus on a lim- 
ited phase space, which allows us to explore the complex 
dynamical behavior as a function of system parameters 
and mass loss rates. We consider only point masses, and 
limit ourselves to circular binaries hosting initially circu- 
lar, coplanar planets, and a single primary mass and evo- 
lutionary pathway. Initial pla net eccentricity is unl ikely 
to be important, as shown by Mudryk & Wu (2006). 

We use a f ourth order Herm ite integrator based on the 
algorithm of Hut et al. (1995), with an additional mass 
loss term. Mass loss rates and evo lving stellar radii ar e 
calculated from the SSE models of |Hurley et al.| ( |2000| . 
We use the default values for mass loss rates (assum- 
ing solar metallicity) provided in the publicly available 
version of the code. Each run consists of two massive 
bodies, and 100 test particles representing planets. The 
time step for all particles is set to 10 -2 — 10 -3 of the en- 
counter timescale for the closest pair of particles in the 
simulation, where the smaller step is required for good 
energy conservation for the fastest mass loss rate^] The 
stellar mass is updated on the same time step as the 
test particles, and all particles share the same time step. 
When mass loss is included, energy is not conserved by 
default, however in test runs with no mass loss, the code 
conserves energy to 1 part in 10 12 over a 5 Myr run. We 
also confirm that particles on both circular and chaotic 
orbits (e.g. bouncing between the sta rs) conserve their 
Jacob i constants over a 5 Myr run (see Moeckel & Veras 
2012 1. We check that the evolution of the stars and test 



particles follows the expected orbital expansion ( Equa- 
tion]^ up until the point of instability. Note that for the 
realistic (fast) mass loss case, a small eccentricity is intro- 

1 The outcomes remain the same under an order of magnitude 
reduction in the time step. 
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Fig. 2. — Evolution of the the zero- velocity curves for a planet as the primary (cyan) loses mass. The white regions indicate where the 
planet may orbit. The bottleneck surrounding the first Lagrange point opens midway through the mass loss phase, allowing the planet to 
orbit both components of the binary. The dashed-tail shows an episode of orbital bouncing, which continues while the curves remain open. 
By the end of the mass loss phase, the zero-velocity surface has been pinched off, and the planet is trapped around the secondary. A movie 
of this evolution is available at http://www.cfa.harvard.edu/~kkratter/BinaryPlanets/zvelrnovie.rnpg. 



duced into the stellar binary this results in oscillations 
in the Jacobi constant for some surviving test particles. 

4.1. Simulation Details 

We simulate evolving binaries containing primaries 
which evolve from M = 2 — > 0.55M©. Such A stars host 



planets (Johnson et al. 2007), undergo extensive mass 
loss in a Hubble time and are fre quently found in multiple 
systems ( |Raghavan et al.||2010 |. We consider secondary 
masses ranging from mi = U.5— 1.7Mq, and initial binary 
separations from a*,j = 75 — 105 AU. We fix the planet 
separation at 15 AU: far enough to avoid any interaction 
with the primary's envelope, and close enough to be well 
inside the (initial) stability region around the host star. 
Test particles are initialized with zero eccentricity, ran- 
dom phase, on coplanar orbits. In order to expedite the 
calculations, we evolve the planetary and stellar orbits 
with the analytic formulae for the first ~ 1490 Myr, un- 
til the star reaches the beginning of the asymptotic giant 
branch (AGB), at which point the mass loss rate rapidly 
increases (shown in Figure [I]) . 

We consider only the evolution of the more massive 
star, which leaves the main sequence first (the planet 
host). In Section [5] we discuss possible outcomes when 
the second star evolves off of the main sequence. 

In total we explore 70 different binary configurations. 
For each binary configuration, we also explore the influ- 
ence of mass loss rate and collisional cross section. Even 



varying only two parameters produces a complex set of 
outcomes because the probabilistic fate of a planet de- 
pends not only on the initial and final configurations, but 
also on the rate at which they pass through each state, 
which is determined by the mass loss rate and expansion 
of the primary's atmosphere. We run our fiducial cal- 
culations for 30 Myr following the end of the mass loss 
phase, or 10 orbits, as most planets that survive to 10 
M yr remain out to 100 M yr i n several test cases. As seen 
by |Rabl fc Dvorak| ( 11988) and |Holman fe Wiegert| ( |1999| ), 
in the case with no mass loss, most of the orbital evolu- 
tion occurs on timescales of less than 300 bin ary o rbits. 
We discuss these timescales further in Section [5T3l 

Even though all the particles begin with nearly identi- 
cal Cj (small variations are introduced due to variations 
in the planets' initial phases), mass loss induces a spread 
in Cj because they undergo different chaotic evolution- 
ary paths while mass loss is ongoing. This spread pro- 
duces a range of outcomes for planets in a given system. 

During the course of the simulation we keep track of 
each particles closest approach to either star. Particles 
which penetrate the specified stellar radii are removed 
from the calculation and identified as collisions. For the 
primary, the collisional radius is set to be t he maximum 
of eith er the stellar radius calculated from |Hurley et aL 
(2000) or the star's approximate tidal radius as a WD: 
3.8 x 10 -3 AU. The secondary's radius is fixed at li?© for 
all masses to simplify the interpretation of probabilities. 
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Fig. 3. — Time evolution of the distance and Jacobi constant for 
the same system shown in Figure [2] where one planet has become 
trapped in orbit about the secondary. The vertical line shows when 
the particle crosses the stability threshold defined in Equation 111. 
The thick black line indicates the stellar separation, the purple 
line indicates the planet separation from the primary, and the pink 
line the distance from the secondary. The thick dashed red and 
grey lines indicate the stability boundaries around the primary 
and secondary, R p and R s respectively. The thin black line shows 
the initial analytic evolution predicted for the planetary orbit. The 
top panel shows the particle's Cj (green), the value of Cj at the 
Li point (black), and the analytic value for Cj of the particle 
on a Keplerian orbit (blue dashed), with separation specified by 
Equation |2jl as discussed in Section [3] 

Because we neglect tidal effects, we also identify any 
particle that has suffered a close approach within 3 tidal 
radii of either star (but not a direct collision) as having 
had a tidal interaction. We define the tidal radius to be: 



1/3 



(6) 



Particles that do not suffer tidal interactions or colli- 
sions are identified as bound to either the primary or 
secondary, although typically no planets remain around 
the primarjj^J 

5. RESULTS 

We begin by showing the probabilities of capture, colli- 
sion, and tidal encounters for the fiducial runs with real- 
istic mass loss rates and stellar radii in Figure |4| We then 
describe in turn the importance of the mass loss rate, and 
physical radius in controlling captures and collisions. 

5.1. Outcomes from fiducial few-body simulations 

Figure[4]shows the probabilities for captures, collisions, 
and tidal interactions for 60 of our 70 binary systems. 
The tabulated results for all systems are shown in Table 

As shown in Table[l] there are five systems which retain planets 
around the primary: in one case (77x2 = O.SMqjO,^ = 105AU), 
the system evolves up to the edge of the stability threshold, and 
because of phase variations only some planets become unstable. We 
repeated this run several times and found similar results. For the 
remaining systems, those classified as bound to the primary are on 
bouncing orbits that we expect will ultimately result in collisions: 
they happen to be in orbit about the primary at the conclusion of 
the simulation. 



[TJ in the Appendix. In the figure, each panel shows the 
outcome for a fixed initial binary separation for a range 
of masses. For all systems that become unstable, the 
most likely outcome for a planet is collision with one of 
the stars. 

Several trends are readily observed. Collisions typi- 
cally increase towards higher masses at fixed separation, 
while tidal encounters decrease. Both trends are the re- 
sult of the instability threshold being reached at an ear- 
lier evolutionary state, when the primary has an inflated 
radius; wider encounters become collisions, and encoun- 
ters within a few stellar radii are prohibited. Note that in 
order to collide with either star, planets must be set onto 
low angular momentum orbits. As a result, collisions in- 
volve at least one passage back and forth between the 
stars. 

At larger separations and stellar masses, planets are 
frequently (10%) captured by the secondary star. Cap- 
tured particles come in two varieties: those whose en- 
ergies as measured via Cj dictate that they remain in 
orbit only about the secondary, and those who in prin- 
cipal may pass back into orbit around the primary, but 
show no evidence of such movement over the course of 
the simulation. _These two are distinguished in Figure [4] 
by the sign of Cj — Cj(L{), where only those with posi- 
tive values are truly trapped (closed zero velocity curves). 
Clearly, some with open curves (negative values) do sur- 
vive for the duration of the simulation. This result is in 



line with the findings of Heppenheimer & Porco |l977 ) . 
They attribute this behavior to planets entering into spe- 
cial orbits that avoid the conjunction of apocenter with 
the L\ point. In general, more planets are captured in 
systems whose zero-velocity curves close at the end of 
mass loss. 

Both classes of captured planets typically have pro- 
grade, eccentric orbits (e > 0.4), and rapid, retrograde 
precession in the frame of the binary. Using the method 
of the surface of section, we identify many of these plan- 
ets as having qua si-periodic, though formall y chaotic or- 
bital parameters ( Murray fc Dermott| |2000[). Their long 
term (> lOOMyr) survival remains uncertain. Neverthe 
less, weaker tidal effects may also dissi pate energy, caus 



ing su ch orbits to shrink and circularize ( Moeckel & Veras 



2012) 



The captured planets typically have orbits which graze 
the stability boundary, as shown in Figure [3] Such plan- 
ets are extremely likely to go unstable when the second 
star in the system evolves to become a WD. A secondary 
instability cycle might re-populate the original host with 
planets on eccentric orbits. 

Close approaches, q, within li?* < q < 3i?*, are also 
common ( 10%), we discuss our treatment of these plan- 
ets, and the plausible results of close encounters and 
physical collisions below. 

5.2. Tidal encounters 
By excluding planets with previous close approaches 

(q ^5 3Rq) from our calculated probabilities for both cap- 
ture and collision, we expect that our neglect of tidal ef- 
fects is unlikely to bias our results. When particle orbits 
are chaotic, they do not evolve smoothly to ever closer 
values of q, but rather make order of magnitude jumps 
in separation from one close approach to the next. To 



G 



10 L 



to 10 1 

O 



10" 
10 c 



fD 10 1 
.Q 
O 



10" 



All Collision"; 


- - - Pnlli^innw/m 

V^UIIIJIUI 1 VV/ //f| 




+ + ^tid 


<g<3i? tid ♦— ♦ 


C'j-Cj(ii) 

F 1 1 1 — — 1 J — [ 








• / 

*^ + 






y * * 






• • 


+ + 


• 


=80.- 

1 1 1 1 ' 1 w 


a *,»= 85 ■ 

! 1 1 1 12 1 





1 1 1 1 [ 


1 1 1 1 1 


1 1 1 1 1 




• 


+ 

+ 


• • 

• 


+ 

a *.* =95 +- 

i i i ■ ™ i X- 


. . . . *a,,=iqp 




+ 

a* ^=105 


0.8 1.0 1.2 1.4 1.6 

m 2 


0.8 1.0 1.2 1.4 1.6 

m 2 


0.8 


1.0 1.2 1.4 1.6 

m 2 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 

1.0 
0.8 
0.6 
0.4 

0.2 ^ 
0.0 
-0.2 
-0.4 



Fig. 4. — Probabilities for collision, capture, and tidal encounters in our fiducial models. Masses are in units of Mq and separations in 
AU. Each panel represents a different semi-major axis for the binary, while different secondary masses are shown along the x-axis of each 
panel. All collisions (black), and collisions with the primary (red dashed) typically increase with increasing mass of the secondary, and 
decrease with increasing binary separation. The decrease in tidal encounters (cyan plusses) with increasing secondary mass is due to more 
collisions with the AGES atmosphere due to the earlier onset of chaotic orbits. Whether or not the zero-velocity curves have pinched off, 
preventing continued planet bouncing, is indicated by the sign of Cj — Cj(Li) (blue diamonds). Positive values show that the curves are 
closed. 

demonstrate this behavior, we show in Figure[5]the evolu- 
tion of the closest approach distance for a small sample 
of particles from one run. We see that a particle that 
reaches a separation of a few stellar radii has likely not 
had previous encounters at comparable distances, and 
thus not experienced significant tidal damping. For plan- 
ets that suffer encounters within a few stellar radii, tidal 
disruption is the most likely outcome. 



5.2.1. Tidal disruption versus capture 

Planet encounters at many stellar radii lose almost no 
energy because the energy dissipated scales as a~ 6 , where 
a p is the periastron. When tidal effects become impor- 
tant at a few tidal radii, close approaches sap energy 
from the planet's orbit. In principal these encounters 
can decrease the planet's velocity (and thus increase Cj), 
enough to bring it from an unstable to stable orbit in a 
binary (or from unbound to bound in a single star sys- 
tem). Tidal capture has been suggested as a formation 
mechanism for X-ray binaries through close encounters 
between stars and compact objects in glob ular clusters 
dFabi an et al.|1975 1, and for close-in planets (Faber et al. 

mbl. — 
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Time (Myr) 
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Fig. 5. — Closest approach q in units of solar radii for a sample of 
test particles in one binary system. Each color represents a differ- 
ent particle trajectory. Note that particles typically undergo large 
jumps in close approach, suggesting that their orbits are chaoti- 
cally sampling the orbital parameter space, not decaying to ever 
smaller separations. 



Recent hydrodynamical simulations suggest that close 
approaches lead to the disruption or ejection of a Jupiter- 
mass planet, rather than long-term survival (Guillochon 



et aLpOlT ). At relatively larger distances, energy is dis- 
sipated through tidal heating of the planet, causing its 
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orbit to shrink. As the orbit decays, the planet begins 
losing mass through tidal stripping; this stripping and 
mass transfer typically lead to the ejection of the planet 
to larger separations at high velocity, or to total plane- 
tary disruption. 

The binary scenario is somewhat different than the 
case of a single-star planet interaction. Tidally stripped 
planets that would have been ejected in the single-star 
scenario could still be bound to the binary system, and 
therefore continue to evolve chaotically following their 
"ejection" in the tidal encounter. In this case, the planet 
might survive or suffer a collision just as any other planet 
in the system. These "shaved" , heated planets might be 
much brighter than their cooling age. A direct physi- 
cal collision with a MS star, (in our case, the secondary) 
should also lead to destruction, as the tidal radius of 
Sun-like stars is comparable to their physical radius for 
Jupiter-like planets. Physical collisions might also be di 
rectly observabl e with upcoming transient surveys ( Met 
zger et al.||2012| ) 



In the case of an encounter with the primary as a WD, 
the tidal radius is much larger than the physical radius of 
the WD; thus disruption is likely prior to a direct physi- 
cal collision. Disruption can lead to the formation of an 
accretion disk around the compact object, and the infall 
of material onto it. Such debris disks around WDs po- 
tenti ally serve as a sourc e for long lasting metal pollution 
(e.g. Klein et al. 2011 and references therein). Similar 
encounters have been studied in the context of the tidal 



disru ption of a star by a massive black hole (e.g. Rees 
19881. 



5.3. Collisions with AGB Atmospheres 

AGB stars spend a short time (~ 10 5 yrs) inflated to 
radii of order 1 AU. Consequently we find that the ra- 
dius distribution of collisions with the primary is highly 
bimodal: they either occur at radii near the maximum, 
or the minimum radius in our models-the tidal radius of 
the WD (see Figure [T]). Fi gure shows the distribution 
of collision radii and times for alTcollisions in our fiducial 
models. We measure time in units of the host binary's 
orbital period at the onset of instability, and t = is set 
to the moment when the system crosses the instability 
threshold. We can identify several trends. First, the bi- 
modality of the distribution is easily observed: only 0.1% 
of collisions occur at intermediate radii. Secondly, we can 
see that closer binaries are destabilized earlier, and there- 
fore collide with the AGB atmosphere at slightly smaller 
sizes while the atmosphere is still expanding (note the 
blue points are shifted to smaller radii compared to the 
green and red points). Lower mass secondaries produce 
AGB collisions only at close separations, whereas more 
massive secondaries drive the system unstable quickly 
enough to drive AGB encounters at all binary separations 
considered. We see that most collisi ons occur within the 
first ~ 200 binary orbits, as seen by Holman fc Wiegert] 
( |1999[ ) and |David et al] ( |2003[ ). We discuss collision 
timescales in more detail in the Section 15.41 

There are several plausible physical outcomes following 
collisions between the planet and AGB star. Collisions 
with a partially evolved star might lead to the engulfment 
of the planet in the stellar envelope, and ultimately to its 
in-spiral and disruption inside the star. In the process, 
the planet might transfer energy, angular momentum, 



and heavy elements into the stellar envelope. 

Collisions with the more loosely bound AGB atmo- 
sphere likely lead to a qualitatively different outcome. 
Such interactions might form a low mass cataclysmic bi- 
nary system as the planet penetrates the AGB envelope 
before its ejection and accretes from it (ILivio fc Soker| 



19841) . It can also affect the structure of the planetary 
nebulae that ultimately forms around the WD remnant 
(Soker 2001), in a similar manner to stellar binary com- 
panions. Trie planet might also completely eject the en- 
velope, depending on the its penetration depth, some- 
what similar to the case of a common envelope scenario. 
In this case, kinetic energy dissipated in the collision and 
envelope ejection could be sufficient to capture the sur- 
viving planet into a close orbit around the stellar rem- 
nant. This latter scenario provides a unique possibility 
for placi ng planets in tig ht orbits around white dwarfs 
(see also|Bear et al. 2011 for a related discussion). Note 



that if tidai dissipation leads to orbits near the tidal ra- 
dius of the W D, such syst ems would reside in the WD 
habitable zone ( |Agol||2011[ ). 



5.4. Varying mass loss rates 

Mass loss due to stellar winds occurs on a timescale 
that is slow (adiabatic) compared to the planet and stel- 
lar orbital period. However, there is another relevant 
timescale in the calculation: the average time for a planet 
to reach its final "outcome" after becoming unstable, i.e. 
collision or a quasi-stable orbit. If the stellar mass and 
binary separation evolve on this timescale, then the prob- 
ability for a collision or capture is determined by the in- 
tegrated probabilities at each state over the time spent 
at each state. 

In order to isolate the effect of system mass ratio and 
separation, and thus determine probabilities as a func- 
tion of system parameters exclusively, we run two sets 
of comparison simulations with unrealistically long mass 
loss timescales set to a constant mass loss rate over 2 
Myr, and 20 Myr. We fix the collision radius to the 
solar radius at all times. In both cases, we compute 
the timescale over which 75% of particles destined to 
collide suffer a collision; this serves as our "equilibrium 

timescale," and ranges from 0.2 < t co u < 1 Myr, depend- 
ing on the system. Changing the fraction of collisions 
which defines our equilibrium timescale by 10% has min- 
imal effect on the timescales. Thus only the 20 Myr mass 
loss timescale case evolves adiabatically with respect to 
all timescales in the problem. 

In this limit, nearly all (99%) of particles suffer a col- 
lision in all of the systems that become unstable. That 
collisions increase in the long mass loss timescale limit 
is easily seen from the evolution of the planets' Cj. At 
the onset of instability, the L\ bottleneck is near its peak 
width. This facilitates exchanges between stars, which in 
turn are responsible for setting planets on collisional or- 
bits. Capture occurs as the bottleneck closes off, which 
only happens at the tail end of mass loss. In the case 
where the mass loss timescale is long compared to the 
collisional time, all of the particles will undergo suffi- 
cient exchanges to collide prior to the closing of the zero- 
velocity curves, and therefore few are left at the end of 
the evolution. We now use these long timescale runs to 
show that the relative collision probability between the 
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Fig. 6. — Distribution of collision times and radii for the fiducial set of runs with realistic radii and mass loss rates. The inset panel is 
zoomed in on the radius range of the AGB atmosphere. Masses are in unites of Mq, and separations in AU. Color indicates the initial 
binary separation, and symbol size the mass of the secondary. Time is measured in units of each host binary's orbital period, with t = 
set to be the point at which the system violates the instability criterion in Equation (TlJ) . Collisions with the primary arc indicated by 
circles, while those with the secondary are shown with squares. A movie which shows tne population of this phase space is available at 
http://www. cfa.harvard.edu/~kkratter/BinaryPlanets/radcollmovie.mpg. 

primary and secondary scales inversely with mass ratio. 
5.5. Collision Rates: Long mass loss times 

Once planets become unstable, in the process of bounc- 
ing back and forth, they sample the orbital parameter 
space about each component, defined by the area en- 
closed within the zero velocity curves. We find that the 
relative collision rates between the stars are well corre- 
lated with the area of the stable region around each star 
as defined by the critical radius in Equation . When 
the stability radius R s or R p is smaller, the physical stel- 
lar radius subtends a larger fraction of the allowed orbital 
area, and thus there is a higher probability for suffering 
a collision. The area ratio depends on the square of the 
stellar mass ratios. Collision probabilities scale inversely 
with the area around the star as shown in Figure[7j mean- 
ing that collisions occur predominantly with the star that 
is less massive. 

In these runs, the stellar radii are fixed at IRq 
throughout the calculation. We also run the same models 
with the collision radius reduced by a factor of 5, and find 
that the collision probabilities remain at 99% across all 
unstable systems, although the collision times are shifted 
to slightly larger values, as shown in Figure[8j Our inter- 
pretation of this trend is that particles set on collisional 
trajectories have nearly radial orbits. This property, and 
the corresponding high collision rate might well be the 
result of considering only coplanar planets. 



o 
U 
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0.45 0.50 0.55 0.60 0.65 0.70 0.75 

R s /( R p + R s ) 

Fig. 7. — Collisional probability for m\ and 7712 as a function 
of the relative area of the secondary stability region R s to the to- 
tal available stable area: R^/(Rp + Rg). These probabilities are 
calculated for the simulation set with slow mass loss (20 Myrs), 
and constant radius. The lower bound on area for which we show 
collisions illustrates the minimum relative areas achieved for com- 
parable primary and secondary masses. 



5.6. Free floating planets: ejection 

In the results presented here, destabilized planets ei- 
ther collided with one of the binary components, or ex- 
changed hosts. No planets were ejected from the system. 



This is because the minimum Cj reached in the systems 
always prohibited escape. For lower Cj, zero- velocity 
curves take on a "c" shape, opening up to the circumbi- 
nary orbital region, thus allowing escape. In several ex- 
ploratory integrations of non-coplanar systems (i ~ 10°) 
and systems with more unequal initial binary mass ra- 
tios (/1 ~ 0.1), we found that planetary ejections are 
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Fig. 8. — Distribution of collision times for the slow mass loss 
(20 Myrs) runs, with collision radius set to Rq (red), and Rq/5. 
(green). As expected, the collision rate peaks slightly earlier for 
larger stellar cross section. Nevertheless, the relative similarity 
demonstrates that colliding particles are on sufficiently low angular 
momentum orbits that collisions destined to reach a stellar radius 
will penetrate to l/5th of that distance. 



possible. For lower mass secondaries, lower Cj may be 
reached because the secondary and planet can be more 
tightly packed prior to instability. Non-coplanar orbits 
likely allow for lower Cj through an expanded phase 
space in which eccentricity (and inclination) excitation 
increase the velocity of the test particle at a fixed loca- 
tion. We defer further discussion of escaping planets to 
future work. 

6. FREQUENCY OF DESTABILIZED PLANETS 

We have thus far shown that the planetary TEDI can 
produce a wide range of astrophysically interesting con- 
figurations. We now try to estimate the occurrence rate 
of planet destabilization. At present, there are few con- 
straints on both the overall frequency of planetary sys- 
tems in binaries, and on wide orbit systems in general. 
Nevertheless, we can calculate the frequency of desta- 
bilized systems for various planetary orbital distribu- 
tions. Given the simplified approach we take, the results 
should be interpreted with caution; they are dependent 
on highly u ncertain model assumpti ons. 

Following Perets & Kratter (2012 1 we construct a pop- 
ulation of stellar binaries, whose properties are chosen 
randomly fro m the appropriate obse rved distributions 



for A-G stars (Raghavan et al. 2010). The periods are 
chosen from a distribution of orbital separations that 
is log-normal. The mass ratio of the binary, TO1/TO2, 
is chosen from a Gaussian dist ribution (mean of 0.6 
and d ispersion of 0.1). Following Fabrycky & Tremaine 
(2007), the dist ribution of initial eccentriciti es is a func- 
tion of period (Duquennoy & Mayor 1991). For peri- 
ods shorter than 1000 days, the eccentricity was chosen 
from a Rayleigh distribution [dp cx eexp(— (3e 2 )de\ with 



(e 2 > 



1/2 



/3-V2 = o.33. For periods longer than 1000 



days, the eccentricity was chosen from an Ambartsum- 
ian distribution (dp — 2ede), corresponding to a uniform 
distribution on the energy surface in phase space. 
We populate the binary systems with single planets. 



Since the distribution of planet orbital properties is un- 
known, we consider three cases : (1) period distributions 
similar to that of low mass binaries (log-normal), (2) 
period distributions similar to that of high mass bina- 
ries (log- uniform), and (3) uniform distribution of peri- 
ods. T he stability criteria derived by |Holman fe Wiegert| 
(1999) does not account for the planetary eccentricity, 
and we therefore do not specify it. 

For the purposes of calculating frequencies, we inte- 
grate our planet distribution down to 3Rq from the pri- 
mary star, however we consider all planets at initial sep- 
arations below 7 AU to be "stable", in that they may 
interact with the stellar envelope or wind, and thus will 
not follow t he simple adiabatic mas s loss evolution envi- 
sioned here ( Villaver fc Livio|[2009 ) 



For all w i der sy stems, we use Equation (1) in Holman 
fc Wiegert| ( |1999[ ) (which includes the binary eccentric- 
ity dependence) to verify whether the initial planetary 
system is stable at birth, and if it is not, we generate 
a replacement system. We construct 1000 initially sta- 
ble planetary systems for each primary mass bin ranging 
from 1 — 3M for each of the three planetary orbital 
distributions. Each primary mass bin is 0.1 Mq. 
We then use the initial to final mass relation deter- 



mined by Salaris et al. (2009), to find the final mass of 
the primary star as it becomes a WD. Using Equation ([2| 
to calculate the evolved orbits of the binary and plane- 
tary companions, we check whether the planetary system 
has crossed the instability threshold. 

The total fraction of unstable systems is recorded; the 
results are shown in Fig. [9j This figure provides an es- 
timate for the fraction of destabilized planetary systems 
for a given primary mass host. Only primaries more mas- 
sive than IMq become WDs in a Hubble time, and we 
therefore use this as a lower limit. Overall we find that 
a few up to ~ 10% of all WD binary systems that con- 
tained planets at birth might undergo the TEDI. 
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Fig. 9. — The frequency of destabilized planetary systems as a 
function of the primary host mass. Three cases are shown, corre- 
sponding to a log-normal, log-uniform and uniform distribution of 
planetary orbital periods. 



7. SUMMARY 

In this paper we have examined the effect of binary 
stellar evolution on the fate of wide-orbit planets. We 
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have shown that a substantial fraction of binary plane 
tary systems may be susceptible to the planeta ry version 



of the Triple Evolution Dynamical Instability (Perets & 
Kratter ||20l2l ) , where mass loss driven orbital expansion 



causes planetary orbits to become unstable and chaotic. 
We have considered circular binary systems and coplanar 
planets so that we may exploit the simplicity provided by 
the circular restricted three-body problem. 

As a result of the chaotic orbits induced by mass loss, 
planets can suffer close encounters or collide with ei- 
ther component in the binary. When close encounters 
occur while the mass losing star is on the asymptotic 
giant branch, dissipation of energy in the envelope can 
potentially lead to planetary captures on close orbits. 
Such planets might reside near the tidal radius of the 
white dwarf following envelope ejection; these planets 
would fall within the white dwarf habitable zone, and 
are potentia lly discoverable with future transit surveys 
(Agol 2011). We note that tidal heating might provide 
a source of geother mal energy to such p lanets. Recent 
numerical work by Kim & Taam (2012) also finds that 
the mere proximity of a planet to the low density atmo- 
sphere might excite spiral arm structure observable with 
ALMA. Encounters within a few tidal radii of either the 
main-sequence secondary, or the newborn white dwarf 
likely result in planetary disruption. Planetary disrup- 
tion around the white dwarf might contribute to the ob- 
served atmospheric pollutio n ( Farihi ct aL 20 09] [Dufour l 
et al.|2010| [Klein et al.|2011||Bochkarev fc Rafikov|201l} . 

Planets that do not sutler close stellar encounters can 
become captured by the secondary. Mass loss and planet 
hopping drive a fall and then rise in the planets' Jacobi 
constants. Mass loss first allows planets to hop back and 
forth between the stars; in some cases the late time rise 
in the Jacobi constant causes planets to become captured 



in orbit about the secondary on quasi-stable orbits (see 
Figure[2}. Because such planets reside near the boundary 
of instability, when the lower mass star in the binary 
evolves, a second round of instability is possible. 

We suggest that future surveys for planets in evolved 
binary systems might find: 

• planets around the main-sequence star in a WD- 
main-sequence binary at the edge of orbital stabil- 
ity 

• planets around the WD in a WD-main-sequence 
binary on sufficiently close orbits that they would 
have been engulfed in AGB atmospheres in stan- 
dard evolutionary models 

• planets around the more evolved WD in a WD- WD 
binary at the edge of orbital stability. 

Future work must constrain in detail the outcomes of 
collisions between planets and AGB atmospheres to make 
more detailed predictions for the types of orbits produced 
by this interaction. Similarly, an exploration of multi- 
planet, and non-coplanar systems is required to provide 
more robust observational predictions. 
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and Computation fellowship through the Harvard Col- 
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APPENDIX 

We include in tabular form the outcomes for all binaries in our fiducial models. We classify particles according to 
three possible outcomes: bound, colliding, and tidally interacting with one of the two components. Note that in several 
cases, the total number of interactions sums to greater than the number of particles in the simulation (100). This is 
due to the fact that several planets suffer tidal encounters with both components. We discuss these results in Section 

o 



12 



TABLE 1 

Tabulated results from the fiducial suite of integrations with 100 test 

particles in each binary. particles are classified as either bound, 
colliding, or tidally interacting with either the primary or secondary. 



Mass (M©) Sep. (AU) Bound 1 Bound 2 Coll 1 Coll 2 Tid 1 Tid 2 
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